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Abstract 

A numerical study of the properties of Gaussian pulses propagating in 
planar waveguide under the combined effect of positive Kerr-type nonlinear- 
ity, diffraction in planar waveguides and anomalous or normal dispersion, is 
presented. It is demonstrated how the relative strength of dispersion and 
diffraction, the strength of nonlinearity and the initial spatial and temporal 
pulse chirps effect on the parameters of pulse compression, such as the maxi- 
mal compression factor and the distance to the point of maximal compression. 
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1 Introduction 



A compression of optical pulses in Kerr-type nonlinear media have been subject 
to investigation for many years and continues to attract a certain attention Jj], 
In single-mode fibers with anomalous group- velocity dispersion (GVD) and positive 
nonlinearity the pulse compression is based on the mechanism of higher-order soliton 
generation ||. In single-mode fibers with normal GVD the pulse compression can 
be obtained in the configuration with a grating pair @, [5]. In both cases the self- 
phase modulation (SPM) induced by an intense pulse is used. However, the intense 
pump pulse propagating together with a weak probe pulse can also cause pulse 
compression by the mechanisms of the so-called cross-phase modulation (XPM) ||, 
or the induced-phase modulation (IPM) |7| . 

A possibility of pulse compression in non- dispersive nonlinear bulk media due to 



another nonlinear effect, that of self- focusing, is discussed in ||-[[nj with the aid of 
the paraxial ray approximation, || pLOfl , and by means of the variational analysis, 
[^]. Still another pulse compression technique that uses the self confinement of 
two-dimensional spatial bright solitons propagating in non-dispersive bulk media is 
mentioned in ||11|| , where the two-beam interference technique is used in order to 



ensure that a filamentation (a splitting of the beam into many sub-beams) does not 
occur. 

Moreover, a simultaneous space-time collapse, which can occur in bulk media 
and in planar waveguides under the combined effect of nonlinearity, diffraction and 
anomalous dispersion, may also be useful for pulse compression [T2|, [T3j . This kind 
of collapse gives rise to short pulses with extremely high optical field [14], [L5| [l6j . It 
is realizable both in the case when dispersion and diffraction have comparable effect 
on pulse propagation and in the more general case when one of the effects above is 
dominating (see fL7|). 

On the other side, the interplay of normal dispersion and positive nonlinearity 
causes quite different behavior of the pulse. In optical fibers where diffraction terms 
are not included it leads to a monotonic pulse spreading. However, the inclusion of 
the diffraction term, which is necessary for a planar waveguide, can lead to a pulse 
compression, as it was described in [T7|, [LB]. Besides, in planar waveguide, normal 
dispersion slows the self-focusing of the pulse and causes a splitting of the pulse into 
two pulses [ 118 , I5 |. The effect of splitting of a pulse was observed also in the bulk 
media |2f. 

In this paper a compression of a pulse propagating in planar, self-focusing nonlin- 
ear planar waveguide in the regime of anomalous and normal dispersion is considered. 
The structure of the paper is as follows. In Section 2, the nonlinear Schrodinger equa- 
tion describing dispersive pulse propagation in nonlinear planar waveguides and the 
parameters of pulse compression are introduced. In Section 3, an estimation of the 
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condition of pulse collapse is made with the aid of the so-called method of moments 
2(J . Numerical results describing the influence of the magnitude of nonlinearity, the 



relative strength of dispersion and diffraction and the spatial and temporal chirp of 
the initial Gaussian pulse on the pulse compression parameters are discussed in 
Section 4. 



2 Basic equations 

It is well known that starting from the Maxwell equations for the envelope U (x, y, z, t) 
of the electric field 

E(x,y,z,t) = U(x,y,z,t)e- i ^ t - n ^ 

propagating along the z axis in a planar waveguide with positive, instantaneous Kerr- 
type nonlinearity, one obtains the 2-dimensional nonlinear Schrodinger equation 
(NSE) : 



if the paraxial and the slowly varying envelope approximations are applied and 
the term V ■ (V-E 1 ), the shock term [[H] proportional to 9 ^ E Q t E ^ and higher-order 
dispersion effects can be neglected. 

In Equation 1, £ = j- is the normalized longitudinal spatial coordinate, £ = ^ 



is the normalized transverse spatial coordinate, r = j^ lZ is the normalized local 



t-f3iz 

CU lldilOVClDC OJJClljlCLl K^KJKJL U.l±±Cllj^, / 

time, a = ^f|^ represents the relative strength of dispersion and diffraction, N = 
PoUoWo^/non2 parameterizes the strength of nonlinearity, (3q = - is a wave number, 
f3 n = are dispersion terms, Zf — f3 n wl is the Fresnel diffraction length, w 
is the spatial width of the input pulse, t is the temporal width of the input pulse 
(i.e., duration of the input pulse), Uq is the peak amplitude of the input pulse, and 
n = no + n2\U\ 2 is the refraction index for the Kerr type nonlinear media. Recall 
that a > corresponds to the normal dispersion and a < corresponds to the 
anomalous dispersion. 

As the initial condition we take the Gaussian chirped pulse which is given by (cf. 
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g 2 (l+iCg) 



U(£, r, C = 0) = e —e * , (2) 

where C% (C T ) is the spatial (temporal) pulse chirp (the focusing spatial chirp cor- 
responds to < and the focusing temporal chirp corresponds to sgn(—aC T ) 
< 0). 
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We will characterize a pulse by its spatial width, ?£>f(C), and the temporal width, 
w t(C), which are defined by 

U(wt,0,() = -U(0,0,C) and U(0, w T , Q = -U(0, 0, C). 

e e 

We also introduce the maximal compression factor 

- T ° 

Cmax / > \ ) 

mm J 

where u> rmm (Cm) is the minimal temporal width of the pulse (see [[18], |23fl). In 
following we call ( m the position of the minimal pulse width. 

Solution of NSE (Equation 1) with the initial condition given by Equation 2 can 
describe a propagation of a dispersive Gaussian pulse in nonlinear planar waveguides. 
It is worth remarking that for the anomalous dispersion regime a solutions of this 



equation can also describe a dispersion-less elliptic Gaussian beam, |24|, £5[ (i.e. a 
cw beam with elliptic Gaussian transverse profile) propagating in a nonlinear bulk 
media. 

In this paper we refer to the case of a = —1 as the cylindrically symmetric 
spatiotemporal pulse; the case of a ^ —1 is to be referred to as the asymmetric 
spatiotemporal pulse. 

In the particular case of the cylindrical spatiotemporal pulse a simple analytic 
solution of the NSE exists which describes a behavior of beam propagating in non- 
linear media by means of the variational approximation [JTJJ] or by means of the 
scaled complex rays formulation within the so-called ABCD matrix formalism (see 
[ pfi| , |27|). For the asymmetric spatiotemporal pulse only a semi-analytical approach 



of p4| , p5[ is known in the literature. 

It is known that some solutions of the two- or three-dimensional NSE can de- 
velop into a singularity of the electric field when the initial pulse power exceeds a 
certain critical value |14| . This phenomenon, known as a pulse collapse, can occur 



simultaneously in space and time for a pulse propagating in planar waveguide with 



the anomalous GVD |L4|, [15|], and also for a dispersion-less beam propagating in 
self-focusing bulk medium. This singularity, however, is obviously non-physical, for 
it emerges just as an artifact of the paraxial approximation made when deriving 
the NSE. In order to avoid this limitation, either a non-paraxial treatment of the 



process of self-focusing or some other effects, such as the nonlinear absorption 
and the saturation of the nonlinear refractive index, should be taken into considera- 
tion. From another hand, the appearance of a non-physical singularity in numerical 
simulations based on NSE can serve as an indication to the real collapse taking place 
in the certain point of space. This is in fact the criterion used in Section 4. 

Studying the details of developing the pulse collapse we leave beyond the scope 
of this paper. Instead, our task is to determine the values of the parameters a and 
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N 2 for which the pulse collapse can occur. For this purpose the so-called method of 
moments p9 could be used. However, it gives only an estimation of the sufficient 



conditions of the pulse collapse, whereas the latter can occur, in fact, at the earlier 
times or on the shorter propagation distances EDI. More precise conditions will be 
obtained by means of the numerical simulations presented in Section 3, (cf. also 



3 Sufficient conditions of pulse collapse 



The method of moments originates from the paper of Vlasov e.a. |[29|| . It can be 
used as an approach to the determination of whether a given initial wave pulse can 



collapse to a singular point in a finite period of time p0[ . An application of the 
method of moments to the NSE may be found in [pC[| . 

In order to formulate the condition of collapse in terms of the strength of non- 
linearity, iV 2 , and the relative strength of dispersion and diffraction, a, we first 
introduce the second moment of intensity 



oo oo 

/(o = / / (e+f 2 )\u\ 2 ^df, 



— oo — oo 



where U is a solution of the NSE given by Equation 1, with the normalization 

f =(-<7)-4r, (<7^0). 

Parameter / can be interpreted as effective beam size measuring the size of the 
area to which most of the energy is confined. 

Assuming that U decay suitably as r — > oo, one can obtain PJ 



where E is the Hamiltonian of the NSE 



E 



AE 



(3) 



fJlk 




2 1 
+ 2 








df 



--iV 2 |f/(£,f)| 4 j^/r 



Because E remains constant during a pulse propagation, i.e. it is independent of Q, 
Equation 3 may be integrated twice to give : 

J(C)=2£C 2 + /(0)C + /(0), 



where I = %. 
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If the right-hand side of the above equation vanished, then the pulse width (both 
spatial and temporal) will decrease to zero in a finite distance leading to beam 
collapse. Therefore a sufficient condition for collapse can occur if the following 
conditions are satisfied ||20|, |22|, BDI : 



E < 
E = 
E > 



and j(0) < 0, 

and 1(0) < -JbEI(0). 



(4) 



For Gaussian input pulse, given by Equation 2, Hamiltonian E can be expressed 
in the following form 



E 



e 



2-2 \ -f 2 -af z 

(j t ] e e 



d^dr 



1 

4 



-<j7rfl 



TV 



In the particular case of flat phase front, = C r = 0, we obtain that 1(0) = 
and because of this two last criterion in Equation 4 are not satisfied. The first 
criterion, E < 0, yields 

N 2 > 1 - a. (5) 

Equation 5 may be considered as the sufficient condition of the pulse collapse in 
terms of the strength of nonlinearity, N 2 , which is proportional to the peak ampli- 
tude, \Uq\ 2 , and the relative strength of dispersion and diffraction, a. The magnitude 
of the parameter N 2 which is sufficient for the pulse collapse to occur increases lin- 
early with \a\. This is not unexpected because the collapse of the pulse occurs 
when the self-focusing caused by the nonlinearity dominates over the broadening of 
a pulse, which is due to the diffraction and dispersion. It is obvious that for smaller 
values of the parameter |a| is, the influence of the dispersion on the pulse broadening 
is weaker. 

Note that the sufficient conditions of pulse collapse can be formulated also in 



terms of the critical initial power P c of the pulse as follows [22, 25 



Pc 

Po 

where P c (a) = J \U\ 2 dfd£ = n\U \ 



M ' 



> 1, 



Pq = 2tt is the initial power of the cylindri- 



cally symmetric pulse (i.e. a = 1). 

We conclude that the decrease of the parameter \a\ leads to the decrease of the 
critical amplitude, \Uo\ 2 , and to the increase of the critical power, P c . 

Note, that the collapse criteria obtained with the aid of the method of moments 
for the particular case of the spatiotemporal symmetric pulse agrees with the result 
of the variational approximation in JTU], [TJJ. 
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4 Numerical results and discussion 



In this Section, the results of numerical solution of the 2+1 dimensional NSE by 
means of the well-known Split-Step Spectral Method (SSSM) |H| with the two 
dimensional (2d) Fast Fourier Transform [521 are presented. The calculations on 
the two-dimensional grid with 512 x 512 points (transverse steps, A£ = At = 0.08) 
and with the longitudinal step depending on the nonlinearity so that for iV 2 = 1, 
A£ = 0.01, were performed. Because of the lack of spatial-temporal cylindrical 
symmetry of the problem it is not possible to simplify calculations by reducing the 
2d Fast Fourier Transform to the one dimensional Hankel Transform developed in 
[j33[ , [34| . Several checks of our numerical procedure were made, which include a 
simulation of beam propagation in the absence of group- velocity dispersion (a = 0), 
repeated testing with different transverse grid and longitudinal step length, and the 
monitoring of pulse energy during each simulation. The latter was kept constant 
with an error < 0.00005 

As initial conditions in numerical calculations we take Gaussian pulse given by 
Equation 2. First, for the case of anomalous dispersion regime we will compare 
conditions of pulse collapse predicted by method of moments with those obtained 
from numerical calculations. Further, with the aid of numerical calculations we will 
study influence of the strength of nonlinearity, relative strength of dispersion and 
diffraction and spatial and temporal pulse chirps on pulse compression parameters. 
Above analysis will be perform both for anomalous as well as normal dispersion 
regime. 

4.1 The anomalous dispersion regime 

In this section the influence of the parameters a and N 2 on the pulse collapse and 
compression will be considered. 

In Fig. 1, a comparison of the conditions of the collapse of pulse predicted by 
the method of moments with those obtained by numerical calculations is presented. 
In our numerical procedure the occurrence of pulse collapse was identified with the 
discontinuity of the phase 0(0, 0,0) in the central point of the pulse u = |w|e^, 
and with a non-monotonic behavior of the intensity in the central point of the pulse 
after reaching the collapse point. The results of numerical simulation are plotted by 
two kinds of points corresponding to the cases when, respectively, the pulse collapse 
occurs or does not occur. The prediction of the method of moments is given by the 
straight line N 2 = 1 — a, (see Equation 5). The boundary line between the collapse 
and the no collapse regions, obtained from the numerical data is approximately 
described by N 2 pa 0.85 — a. It is parallel to the straight line predicted by the 
method of moments, unless the absolute value of a is too small. 
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Therefore, for both methods, the magnitude of the parameter N 2 which is suf- 
ficient for the pulse collapse to occur increase linearly with \cx\. The discrepancy 
appears due to the theoretical idealization of the picture of the collapse where all 
the energy of the pulse goes to the singularity point. This also explains why condi- 
tions of numerical collapse are typically softer than those predicted by the method 
of moments described in Section 3. 

Studying details of pulse collapse we leave beyond task of this paper. Instead, 
we will study the influence of the relative strength of dispersion and diffraction, 
the nonlinearity and the spatial and temporal chirps on the parameters of pulse 
compression under the condition that the pulse collapse does not occur. 

Fig. 2 and Fig. 3 represent the results of calculations of the influence of the 
relative strength of dispersion and diffraction, a, on the maximal compression factor, 
c max , and on the position of the minimal pulse width, ( m , for different values of the 
strength of nonlinearity, N 2 and for Gaussian initial pulse with flat phase front. 
As it could be expected the parameters of pulse compression, c max and Cm,i ncrease 
monotonically with the increase of iV 2 and the decrease of a until collapse conditions 
are reached. This behavior is obvious from the fact that increase of N 2 cause increase 
of pulse self-focusing, it helps to concentrate pulse energy in the center, in addition, 
decrease of a cause decrease of dispersion broadening of the pulse. 

In Fig. 4 the results of numerical simulations of the influence of the initial spatial, 
C% = C, and two cases of temporal, C T = ±C, chirps on the pulse compression 
parameters are presented. In order to distinguish between the above two cases we 
introduce a parameter 

S = sgn(—C^C T a), 

which equals 1 for the case of focusing (defocusing) temporal and spatial chirps and 
equals -1 for the case of focusing (defocusing) temporal and defocusing (focusing) 
spatial chirps. 

As it could be expected, the focusing spatial and temporal chirps, C < 0, S — 1, 
cause the increase of the pulse compression parameters. The explanation is that 
a defocusing chirp spreads the energy out from the center of the pulse, whereas a 
focusing chirp concentrates it there. As the result, the nonlinearity-induced phase 
curvature of the field is, respectively, reduced or enhanced. Similar effect of the 
focusing chirp of the initial pulse takes place in the region close to the collapse. 
Namely, the focusing spatial chirp can hasten the collapse, whereas a defocusing 
chirp can either delay or eliminate it entirely Wh - 

More interesting is the case of S — — 1 (i.e. the spatial focusing chirp and 
the temporal defocusing chirp occur simultaneously). The increase of the maximal 
compression factor occurs only for the case of focusing temporal chirp, C > 0, 
whereas this is not always true for a spatial focusing chirp C < 0, see Fig. 4. 
One can conclude that the temporal chirp has larger effect on the temporal pulse 
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compression than he spatial one. One can expect the reverse situation in the case 
of spatial compression of the pulse. 



4.2 Normal dispersion regime 

In the case of normal dispersion regime the collapse of the pulse doesn't occur. 
However, due to the the spatiotemporal coupling occurring in nonlinear medium 
when both the diffraction and the dispersion effects take place a pulse compression 
can be obtained [IT 



In this section we will study the influence of the relative strength of dispersion 
and diffraction, the nonlinearity and the spatial and temporal pulse chirps on the 
parameters of pulse compression. 

It is seen from Fig. 5 that the maximal compression factor, c ma x, monotoni- 
cally decreases with a, and increases with N 2 . It is clear because smaller value of 
the parameter a has a weaker influence on the dispersion broadening of the pulse, 
moreover the increase of iV 2 leads to the increase of the spatiotemporal coupling 
and nonlinearity induced phase curvature of the field. In the end both effects lead 
to the temporal compression of the pulse. 

From Fig. 6 it is seen that for sufficiently small values of a the parameter ( m 
decreases with iV 2 . However it appears to be practically independent on a once a 
certain threshold level of iV 2 is reached. This fact was explained in |l7j by means of 



the periodic beam narrowing of higher-order spatial solitons. A different behavior 
takes place at larger values of sigma (a > 0.25). Namely, at first Cm increases with 
N 2 for sufficiently small iV 2 and then it slowly decreases after reaching a maximal 
value at the certain value of N 2 . This behavior is explainable by the fact that at 
small nonlinearities the effects of dispersion prevent a creation of spatial solitons. 

In Fig. 7 the results of numerical calculations of the influence of the initial pulse 
chirp on the parameters of pulse compression c m and ( m are presented. The focusing 
spatial and temporal chirps, C < 0, S — 1, cause the increase of the compression 
parameters (c max ans (max) and this behavior appears to be similar to that which we 
have previously observed in Fig. 4 for the anomalous dispersion regime. However, 
in the case of the anomalous GVD c max grows with C much faster that in the case of 
normal GVD. Namely, for the anomalous GVD the maximal compression factor for 
a chirped initial pulse with C = —2 is three times larger than that for an initial pulse 
with flat phase front (C = 0, i.e. c max (C = —2) = 3 x c max (C = 0). For the normal 
GVD the increase of the c max is rather slow, e.g. c max {C = —2) = 1.1 x c max {C = 0), 
and a saturation of the maximal compression factor occurs for the initial chirps below 
-2 (see Fig. 7). 

Moreover, for the case of C < 0,S = — 1 (i.e. the spatial focusing chirp and 
the temporal defocusing chirp) the maximal compression factor increases only for 
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focusing temporal chirp, whereas this is not always true for a spatial focusing chirp 
Q < 0. 



5 Conclusions 

In this paper, the physical conditions of collapse and compression of dispersive 
Gaussian pulses propagating in waveguide with the positive Kerr-type nonlinearity, 
diffraction and the anomalous or normal dispersion are investigated. 

We determine the values of the relative strength of dispersion and diffraction, 
a, and the strength of nonlinearity, N 2 , for which the pulse collapse can occur. For 



this purpose we first present an estimation given by the method of moments |29 
More precise conditions are obtained by means of the numerical simulations based 
on the (2+l)-dimensional Nonlinear Schrodinger Equation (see Section 4). 

We characterize a pulse compression by two parameters: the maximal compres- 
sion factor, c max , and the distance to the point of the maximal compression, ( m , (see 
Section 2). By means of a numerical simulation we study how these two parame- 
ters depend on the parameters N 2 and a, and the initial spatial and temporal pulse 
chirps. We demonstrate that in the regime of both anomalous and normal dispersion 
the increase of the nonlinearity and the decrease of the relative strength of dispersion 
and diffraction cause the increase of the maximal compression factor. Moreover, in 
the case of anomalous dispersion regime the compression factor is maximal in the 
region of 1 — a ~ iV 2 . 

Furthermore, we observe that the increase of the focusing temporal and spatial 
chirps of the initial pulse lead to the increase of the maximal compression factor, 
Cmax- In the case of the anomalous GVD c max grows with chirp, C, much faster 
that in the case of normal GVD, for which a saturation of the maximal compression 
factor occurs. 

Moreover, the increase of the focusing temporal chirp might lead, even in the 
presence of the defocusing spatial chirp, to the increase of the maximal compression 
factor, c max , whereas the defocusing temporal chirp always leads to the decrease of 
Cmax, even in the case of the focusing spatial chirp. It may be concluded, therefore, 
that the temporal chirp has larger effect on the maximal pulse compression factor 
than the spatial chirp. In reverse, it is expected that the spatial focusing chirp has 
larger impact on the beam with than the temporal chirp, independently on its sign. 
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Figure 1: Comparison of the sufficient conditions for pulse collapse predicted by the 
method of moments (straight line) and numerical calculations (filled circle points 
denote pulse collapse and empty circle points indicate no collapse). It was done for 
initial Gaussian pulse with flat phase front, propagating in a medium described by 
two parameters: the strength of the nonlinearity, N 2 and the relative strength of 
dispersion and diffraction, a. 
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Figure 2: The maximum compression ratio, c max as a function of the relative strength 
of dispersion and diffraction, a < 0, for different value of the strength of nonlinearity, 
iV 2 , and for initial Gaussian pulse with flat phase front. 
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Figure 3: The distance to the point of minimal pulse width, ( m , as a function of 
the relative strength of dispersion and diffraction, a < 0, for different value of the 
strength of nonlinearity, iV 2 , and for initial Gaussian pulse with flat phase front. 
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Figure 4: The maximal compression factor, c max , and the distance to the point of 
the minimal pulse width, Cm, as a function of the initial spatial and temporal pulse 
chirps, = C and C T = ±C, respectively. Spatial focusing chirp occurs for 
C < 0, temporal focusing chirp occurs for C T < 0, a = —0.5 and iV 2 = 1.0. 
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Figure 5: The maximal compression factor, c max , as a function of the relative 
strength of dispersion and diffraction, a > 0, for different value of the strength 
of nonlinearity, N 2 , and for initial Gaussian pulse with flat phase front. 
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Figure 6: The distance to the point of the minimal pulse width, Cm, as a function 
of the strength of nonlinearity, N 2 , for different value of the relative strength of 
dispersion and diffraction, a > 0, and for initial Gaussian pulse with flat phase 
front. 
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Figure 7: The maximal compression factor, c max , and the distance to the point of 
the minimal pulse width, Cm, as a function of the initial spatial and temporal pulse 
chirps, = C and C T = ±C, respectively. Spatial focusing chirp occurs for 
< 0, temporal focusing chirp occurs for C T > 0, a = 0.1 and iV 2 = 2.0. 
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